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FOREWORD 


Determination  of  beach  material  characteristics,  and 
their  variance  across  and  along  the  beach,  is  impu*- U«nt  in 
many  coastal  studies  — in  indicating  the  source  and 
direction  of  movement  of  the  littoral  materials  and,  to 
some  extent,  the  probable  reaction  of  the  beach  to  particular 
wave  and  tide  conditions.  Until  very  recently  there  has 
been  no  sound  basis  for  planning  a beach  sampling  operation, 
samples  being  taken  in  various  ways  and  analyzed  by  various 
methods  with  little  if  any  realization  of  the  statistical 
implications  involved.  Studies  have  recently  beo.n  underway 
which  attempt  to  discover  the  best  statistical  model  to  be 
used  in  beach  sampling  (e.g.,  number  and  size  of  samples, 
depth  of  sampling,  lateral  and  longitudinal  spacing  of 
samples),  and  the  following  is  a report  on  one  of  these 
studies. 


This  report  has  been  prepared  by  Dr.  William  C. 
Krumbein  in  pursuance  of  Contract  DA-i49-055~eng~35  with  the 
Beach  Erosion  Board,  which  provides  in  part  for  research  on 
programs  and  procedures  for  sampling  littoral  materials. 

Dr.  Krumbein  is  Professor  of  Geology  at  Northwestern  Univ- 
ersity. 

Views  and  conclusions  stated  in  this  report  are  not 
necessarily  those  of  the  Beach  Erosion  Board. 
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This  report  is  published  under  authority  of  Public  Law 
79th  Congress,  approved  July  31 , 19U5. 
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STATISTICAL  SIGNIFICANCE 
CF 

BEACH  SAMPLING  METHODS 

by 

W.  C.  Krumbein,  Northwestern  University 


INTRODUCTION! 


Characteristics  of  beach  deposits  are  important  in  geological  and 
engineering  investigations  of  shore  ne  processes.  Almost  every  study 
of  beaches  includes  collection  of  samples  for  determination  of  average 
particle  size,  mineral  composition,  and  other  attributes. 

Most  beach  samples  are  collected  with  one  or  more  of  the  followirg 
objectives  in  mind: 

(a.)  To  obtain  an  estimate  of  the  average  coarseness  or 
fineness  of  the  beach  material  for  descriptive  or 
comparative  purposes. 

(b)  To  obtain  information  on  the  characteristics  of 
material  required  for  beach  nourishment  or  for  the 
development  of  artificial  beaches. 

(c)  To  determine  the  source  rocks  or  source  areas  of 
the  natural  sand  or  gravel  on  the  beach. 

(d)  To  determine  whether  there  have  been  changes  in  the 
sources  of  the  beach  material  within  historic  times 
due  to  man's  modification  of  the  beach  area  or  its 

hinterland. 

Inasmuch  as  the  characteristics  of  beach  material  are  in  part 
dependent  on  physical  processes  operating  in  the  beach  ares,  the 
preceding  specific  objectives  may  be  parts  of  a larger  plan  to  x'elate  the 
observed  properties  of  the  material  to  prevailing  waves  and  currents, 
as  well  as  to  the  distance  that  the  material  was  moved  from  its  original 
source  area. 

Methods  of  sampling  beaches  are  largely  empirical  in  that  few 
studies  have  been  made  to  test  the  effectiveness  of  different  sampling 
plans.  The  purpose  of  this  report  is  to  review  beach  sampling  methods 
in  terms  of  some  statistical  principles,  and  to  suggest  tentative  plans 
for  improving  the  representativeness  in  the  samples.  In  order  to  be  self- 
contained  the  report  includes  some  discussion  of  the  statistical  methods 
used.  The  results  indicate  that  no  radical  revisions  of  current  sampling 


procedures  are  needed,  but  that  recognition  of  explicit  design  elements 
in  the  sampling  plan  can  result  in  greater  reliability  of  the  data  for 

^nmnsrshle  eroendit.ure  of  time  and  pffnrt.. 

- - * — -- 

The  writer  is  indebted  to  numerous  individuals  for  suggestions  and 
data  used  in  this  study.  Messrs.  R.  0.  Eaton*  J.  M.  Caldwell,  J.  V. 
Hall,  Jr.,  and  T.  Saville,  Jr.,  of  the  Beach  Erosion  Board  staff  have 
been  especially  cooperative  in  discussing  the  problem  and  in  placing 
beach  information  at  the  writer's  disposal.  Mr.  Eaton  has  made  many 
specific  suggestions,  and  has  read  the  manuscript  critically  for  its 
engineering  implications. 

REVIEW  OF  BEACH  CHARACTERISTICS 


Sampling  programs  to  determine  specific  properties  of  the  material 
being  sampled  are  most  effectively  designed  when  the  general  character- 
istics of  the  material  are  known.  In  a statistical  sense,  each  set  of 
observations  represents  a sample  drawn  from  some  larger  population.  If 
the  population  is  known  to  be  relatively  homogeneous,  a set  of  random 
samples  may  supply  sufficient  data  for  estimating  population  character- 
istics. If  the  population  shows  natural  subdivisions  or  gradients,  the 
sampling  plan  may  be  adjusted  to  take  such  variations  into  account. 

It  is  common  knowledge  that  beach  material  characteristics  change 
along  and  across  beaches  due  to  variations  in  the  "natural  treatments" 
imposed  by  waves  and  currents.  Waves  wash  up  to  the  berm,  but  seldom 
cross  it.  Seaward  of  the  berm  the  sand  is  shifted  about  by  wave  action, 
but  on  the  dry  backshore,  wind  winnowing  may  selectively  remove  grains 
to  the  dunal  belt.  As  a result  of  these  differing  processes,  the 
beach  as  a whole  displays  a zonal  arrangement,  and  is  commonly  divided 
into  a backshore,  foreshore,  and  nearshore  bottom  zone.  The  width  and 
characteristics  of  the  zones  are  functions  in  part  of  the  amounts  of 
beach  material  available,  as  well  as  of  the  larger  environment  of  partially 
enclosed  bay,  open  shore,  or  barrier  system  of  which  the  beach  is  a part. 
Within  this  broader  framework  the  sampling  plan  is  thus  part  of  a larger 
design  that  includes  the  physical  setting  of  the  beach. 

The  problem  of  designing  a comprehensive  sampling  plan  for  such  an 
interrelated  system  is  beyond  the  scope  of  this  preliminary  study,  in- 
stead, attention  will  be  confined  to  straight  stretches  of  open  beach  as 
being  more  representative  of  conditions  commonly  sampled.  Even  here 
the  problem  is  not  simple.  Cyclical  changes  occur  on  the  beach  seasonally 
and  during  storms.  Low  summer  Derms  are  topped  by  fall  and  winter  storms, 
and  new  berms  adjusted  to  prevailing  wave  conditions  are  developed.  Even 
during  a single  season  the  changes  which  occur  during  storms  in  the  neap 
tide  range  are  noticeably  different  from  those  occurring  in  spring  tide 
ranges.  Likewise  there  is  a difference  in  effects  produced  by  storms 
occurring  during  rising  or  falling  tides. 


These  cyclical  changes  effect  beach  widths,  slopes,  and  textures 
of  the  beach  material.  Some  beaches  may  be  composed  of  sand  in  summer 
and  pebbles  in  winter;  in  others  (as  seme  beaches  along  Lake  Michigan) 
the  beach  may  be  composed  of  pebbles  during  quiet  periods  and  of  sand 
during  storms.  So^ie  beaches  are  characterized  by  patches  of  gravel  on  a 
sand  foundation;  the  position  and  size  of  the  gravel  patches  may  vary 
with  storms  and  seasons.  Any  sampling  plan  which  takes  all  these 
possibilities  into  account  must  be  repetitive,  and  must  extend  from  the 
highest  berm  crest  to  a water  depth  sufficient  to  include  areas  of 
significant  amounts  of  sand  movement. 

The  backshore  beyond  the  highest  berm  is  very  seldom  affected  by 
wave  action,  and  then  only  during  severe  storms  or  breakthroughs.  From 
the  sampling  point  o"  view  the  backshore  and  its  belt  of  landward  dunes 
may  be  considered  a > <?!?. lively  unchanging  part  of  the  beach  as  far  as 
hydraulic  effects  are  concerned. 

Inasmuch  as  emphasis  in  the  present  report  is  on  engineering  uses 
of  beach  data,  the  foreshore  and  nearshore  bottom  will  receive  most 
attention  because  they  represent  parts  of  the  beach  most  continuously 
subject  to  hydraulic  forcer,  dome  aspects  of  the  backshore  are  treated, 
however,  to  provide  a basis  for  comparison  with  the  more  active  portions 
of  the  beach. 


BEACH  POPULATIONS 

It  was  mentioned  that  some  knowledge  of  the  populations  from  which 
samples  are  drawn  is  helpful  in  designing  sampling  plans.  There  are 
other  advantages  also,  such  as  in  the  selection  of  s tatisticarl  models 
for  further  analysis  of  tne  data,  and  in  setting  up  confidence  intervals 
for  estimates  of  certain  population  parameters. 

A considerable  body  of  knowledge  about  some  beach  populations  is 
available  from  the  it...,  y geological,  and  engineering  studies  that  have 
been  made.  In  part  tMs  knowledge  is  used  by  implication  rather  than 
e::pjLi.citly.  Some  of  the  material  is  reviewed  here  to  indicate  that  at 
least  six  kinds  of  distribution  functions  are  commonly  encountered  in 
beach  data. 

In  any  small  area  on  a large  beach  the  average  particle  size,  particle 
shape,  mineral  composition,  moisture  content,  beach  firmness,  and  other 
characteristics  are  relatively  homogeneous,  and  closely  spaced  samples 
commonly  show  approximately  normal  distributions  of  these  properties. 

The  distribution  of  average  particle  size  in  closely  spaced  beach 
sane,  samples  tends  to  be  normal,  although  the  distribution  of  grain  sizes 
in  any  single  sample  is  skewed.  The  individual  samples  may  be  "normalized" 
by  a log  transformation,  as  is  commonly  done  graphically  by  plotting 
the  size  data  on  logarithmic  paper. 
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Some  beach  characteristics  have  a directional  or  angular  sense, 
such  as  foreshore  slopes,  orientations  of  grain  axes,  and  other 
attributes.  These  commonly  show  symmetrical  distributions  which  can 
be  approximated  by  the  normal  curve.  In  theory,  however,  these  angular 
distributions  are  circular  normal • 

Some  beach  data  show  highly  skewed  frequency  distributions  which 
cannot  be  normalized  by  a log  transformation.  These  are  Gamma  distributions, 
and  they  occur  when  the  rarer  chemical  constituents  of  sedimentary 
deposits  are  studied,  for  example.  In  some  cases  the  distribution  of 
bedding  thicknesses  follows  a Gamma  distribution. 

The  preceding  kinds  of  distributions  apply  to  continuous  data.  In 
contrast  to  these  are  discrete  distributions,  where  the  variables  takes 
on  only  integral  values.  These  distributions  aie  encountered  when 
mineral  grains  in  sands  are  counted,  such  as  in  the  analysis  of  heavy 
minerals.  The  distributions  are  developed  by  taking  subsamples  of  fixed 
size  (say  25  to  50  grains),  and  counting  the  number  of  specific  occurrences 
per  subsample.  If  the  mineral  counted  is  fairly  abundant,  each  sub- 
sample will  have  one  or  more  grains,  whereas  if  the  mineral  is  rare, 
some  of  the  subsamples  may  not  contain  any  of  the  grains.  These  con- 
ditions give  rise  to  two  kinds  of  distributions.  The  binominal  distribu- 
tion, relatively  symmetrical,  is  found  when  the  minerals  are  abundant, 
and  the  ^oisson  distribution  occurs  for  very  rare  minerals.  Poisson 
distributions  are  commonly  highly  skewed. 

Table  1 summarizes  the  six  kinds  of  distributions,  and  includes 
several  examples  of  each  type.  The  table  serves  as  background  material 
for  later  mention  of  certain  beach  populations.  Figure  1 shows  an 
example  of  each  kind  of  distribution  for  illustrative  purposes. 

BEACH  ZGNES 


It  was  mentioned  that  beaches  may  be  roughly  divided  into  the  back- 
shore,  the  foreshore,  and  the  nearshore  bottom,  as  indicated  in  Figure  2. 
In  a ctatistical  sense  these  constitute  three  or  more  "sampling  strata", 
each  of  which  has  its  own  characteristic  populations.  The  strata  may 
be  considered  as  separate  units,  they  may  be  further  subdivided,  or  they 
may  be  considered  as  parts  of  a "superpopulation"  which  comprises  the 
beech  as  a whole.  In  taary  instances  there  are  continuous  gradations 
between  the  zones,  so  that  the  population  as  a whole  may  have  gradients 
from  one  zone  to  another. 

An  easily  visualized  example  of  populations  is  illustrated  by  the 
moisture  content  of  beach  sand.  In  the  nearshore  bottom  all  the  pores 
among  the  sand  grains  are  filled  with  water  and  the  moisture  content 
is  of  the  ordar  of  35  to  LiO  percent.  On  the  foreshore  the  samples  close 
to  the  water  line  may  have  20  to  25  percent  moisture,  and  as  samples 
are  collected  higher  on  the  foreshore  this  percentage  rapidly  drops  to 
between  5 and  10  percent.  On  the  crest  of  the  berm  there  may  be  only 
1 or  2 percent  moisture,  and  the  somewhat  lower  elevation  of  the  backshore 
may  develop  slightly  higher  moisture  values. 
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TABLE  1 


FREQUENCY  DISTRIBUTION  OF  BEACH  DATA* 

1 . NCRMA1  . OR  NORMAL  APPROXIMATION 

Average  particle  size 

Particle  roundness  and  sphericity  in  some  samples 
Moisture  content  of  beach  sand 
Beach  firmness 
Porosity  of  beach  sand 

Orientation  of  particle  ares  when  range  is  small 
Beach  slopes  expressed  in  angles 
Percentages  of  abundant  minerals  in  sand 

2.  LOG  NORMAL 

Particle  size  distribution  in  single  samples 
Permeability  of  beach  sands 
Thickness  of  bedding  (?) 

3.  CIRCULAR  NORMAL 

Orientation  of  particle  axes  when  range  is  large 
Orientation  of  cross-bedding  in  beach  deposits  (?) 
Distribution  of  angle  of  approach  of  waves 

iu  BINOMIAL 

Abundant  minerals  in  sand  or  of  rock  types  in  gravel 

(when  data  are  expressed  as  number  of  particles  occurring 
in  subsamples  of  fixed  size.) 

Number  of  frosted  grains  in  some  dune  sands  (in  subsamples  of 
fixed  size). 

5.  POISSON 

Rare  minerals  in  sand  (or  of  rock  types  in  gravel),  when  data 
are  expressed  as  number  of  particles  per  subsample  of  fixed 
size. 

6.  GAMMA  DISTRIBUTIONS 

Thickness  of  bedding  (?) 

Percentage  of  rarer  minerals  in  sand 

Particle  roundness  and  sphericity  in  some  samples 


* Exceptions  to  this  classification  occur.  Note  also  that  percentages 
are  usually  binomial  or  Poisson  basically,  but  as  continuous  variables 
they  apparently  may  be  treated  as  normal  approximations  or  Gamma  dis- 
tributions. 
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Number  Freqtincy  Number  Number  Frequency 


Normal  (Beach  Lrmress) 


Binomial  (Limestone  pebbles  in 


Poisson  (Number  of  magnetite  grains 
in  subs  ample?  of  40 ) 


Circular  Normal  (Orientation  of 


FIGURE  I - DISTRIBUTION  FUNCTIONS  ENCOUNTERED 

IN  BEACH  SAMPLING 
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FIGURE  2 BEACH  ZONES  (SAMPLING  STRATA) 
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In  treating  the  moisture  data,  one  may  consider  that  the  entire 
beach  has  one  superpopulation  of  moisture  content,  with  a sharp  gradient 
from  the  water  line  landward.  On  the  other  hand,  one  may  consider  the 
nearshore  bottom  and  the  back shore  as  two  distinct  populations  with  a 
transition  zone  between.  In  fact,  one  may  set  up  several  parallel  popu- 
Nation  bands  ar.rnsR  t*h6  HpnpnH-j na  upon  the  t, ijjsit-stions  he 

to  impose  upon  any  single  "population." 

Although  the  changes  may  not  be  so  obvious  to  the  eye,  similar 
variations  occur  in  particle  size,  mineral  composition,  beach  firmness, 
and  other  attributes  from  zone  to  zone.  The  degree  of  variability  itself 
varies  with  different  properties,  so  that  the  beach  as  a whole  may  be 
considered  as  a complex  of  populations  of  various  attributes,  each  with 
its  own  characteristics  and  rates  of  change. 

In  part,  the  Sampling  plan  depends  upon  the  population  model  adopted. 
In  many  instances  it  is  preferable  to  treat  each  zone  separately,  although 
in  some  studies  a broad  over-all  value  for  the  beach  ss  a whole  may  be 
adequate. 


CLASSIFICATION  OF  SAMPLING  METHODS 


Methods  of  beach  sampling  have  in  large  part  developed  from  ex- 
perience, and  only  in  realtively  recent  years  has  special  attention  been 
paid  to  the  statistical  problem  of  obtaining  the  best  kind  of  examples 
for  a given  study.  During  the  past  two  decades  statisticians  have 
developed  sampling  theory  to  an  advanced  level  in  many  fields  (such  as 
agricultural  sampling,  opinion  sampling,  etc.).  The  body  of s ampling 
theory  available  applies  equally  well  in  principle  to  other  fields,  and 
some  of  the  definitions  and  concepts  are  cf  direct  value  to  beach  sampling. 

A distinction  should  bs  made  between  purposive  selection  and  random 
sampling.  In  purposive  selection  the  samples  are  restricted  to  deposits 
or  places  considered  by  the  worker  to  be  especially  typical  of  the  con- 
ditions being  ampled.  A random  sample  is  one  in  which  some  elements  of 
random  selection  has  been  introduced  to  avoid  personal  bias,  and  it  is 
in  connection  with  the  latter  kind  of  samples  that  sampling  theory  is 
developed.  Cochran  (1953)  develops  this  point. 

Trie  requirement  that  an  element  of  randomness  be  introduced  into 
sampling  does  not  mean  that  the  sample  is  to  be  taken  blindfolded,  as 
it  were,  without  regard  for  the  characteristics  being  sampled.  Rather, 
the  requirement  is  that  some  process  of  randomization  be  appl ied  to  the 
procedure  of  collecting  the  sample.  This  randomization  process  may  take 
different  forms,  depending  upon  the  kinds  of  samples  being  collected. 

The  distinction  is  brought  out  in  the  following  definitions  of  sampling 
procedures.  The  formal  definitions  are  given  by  Cochran  (1953/;  they  are 
paraphrased  here  to  apply  more  directly  to  beach  sampling,  either  along 
a profile  or  over  an  area. 

In  sampling  sedimentary  deposits,  the  unit  of  sampling  is  a small 
volume  of  the  deposit  r?Jher  than  an  individual  grain,  although  in  some 


/ 
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studies  the  grains  themselves  may  be  sampled.  In  the  following  treat- 
ment it  will  be  assumed  that  the  beach  sample  is  a cylinder  of  sand 
about  3 inches  in  diameter  and  2 inches  deep.  Four  kinds  of  sampling 
procedures  are  available,  which  will  be  illustrated  first  by  samples 
collected  along  a profile  line  across  the  beach. 

1=  Simple  random  sampling  is  a method  of  selecting  a sample  of 
given  size  such  that  every  one  of  the  possible  samples  of  this  size  along 
the  line  has  an  equal  chance  of  being  selected- 

2.  Stratified  sampling  involves  dividing  the  sampling  line  into  L 
non-overlapping  segments  (strata)  which  together  constitute  the  entire 
length  of  the  line.  A set  of  L samples,  representing  one  sample  collected 
at  random  from  each  segment,  is  a set  of  stratified  random  samples. 

3.  Systematic  sampling  is  performed  by  taking  a sample  at  random 
from  an  arbitrarily  defined  segment  of  length  k,  and  thereafter  selecting 
a sample  from  the  same  relative  position  in  all  remaining  segments  of 
length  k along  the  original  sampling  line. 

U.  S+ ratified  systematic  sampling  is  performed  by  first  dividing 
the  line  into  L strata  as  under  (2).  A sample  is  collected  at  random 
in  the  first  stratum,  and  thereafter  a sample  is  taken  from  this  same 
relative  position  in  each  of  the  remaining  strata. 

Figure  3 comprises  plan  views  of  a beach  profile  1,000  feet  long 
extending  from  the  inner  edge  of  the  backshore  to  the  30-foot  depth  in 
the  nearshore  bottom.  The  four  lines  in  the  figure  represent  the  same 
profile  sampled  in  four  different  manners.  In  profile  A five  samples 
are  collected  at  random.  One  way  of  selecting  the  samples  is  to  consider 
the  profile  as  comprising  u,000  consecutive  sampling  positions,  each  3 
inches  long,  numbered  from  0001  to  U,000.  A table  of  random  numbers 
(such  as  Jhose  in  Dixon  and  Massey,  p.  290)  is  entered  at  some  convenient 
point,  and  the  first  five  U-digit  numbers  are  listed,  ignoring  numbers 
larger  than  U,000,  and  omitting  duplicates.  Say  the  numbers  are  1326, 

272li,  02i*9,  1196,  and  3253.  Then  a sample  is  taken  from  the  2ii?th, 

1126th,  1196th,  2?2Uth,  and  3253rd  3-inch  segments,  representing  distances 
of  62.25,  261.50,  299.00,  661.00,  ana  6l3.25  feet  from  the  landward  end 
of  the  profile.  These  constitute  a set  of  simple  random  samples,  with 
some  qualifications  developed  later. 

In  profile  B the  beach  is  first  stratified  into  backshore,  upper  and 
lower  foreshores,  and  say  two  nearshore  bottom  zones.  7hece  are  natural 
strata  suggested- by  past  experience  on  the  beach,  and  they  constitute  the 
L segments  of  the  definition.  Each  stratum  is  now  considered  as  a sampling 
segment,  and  a separate  random  number  is  drawn  for  each  stratum  from  a 
range  of  values  equal  to  the  number  of  possible  3-inch  ssmDle  spacings 
along  the  line.  The  five  resulting  samples  constitute  a set  cf  sti'atified 
random  samples  inasmuch  as  each  stratum  was  sampled  in  one  random  position. 
For  strata  of  different  lengths  it  may  in  some  cases  be  desirable  to  take 
a number  of  samples  from  each  stratum  proportional  to  the  length  of  the 
stratum. 
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FIGURE  3 - SAMPLING  PLANS  (BEACH  PROFILE) 


In  profile  C the  sampling  line  is  considered  as  having  1*,000  possible 
sampling  sites  as  in  A,  but  it  is  now  arbitrarily  decided  to  consider 
the  profile  as  being  composed  of  five  sampling  s egments  of  equal  length. 
These  are  the  k groups,  each  of  which  has  800  of  the  original  h,OOC  3-inch 
sampling  elements.  A single  number  between  0C1  and  800  is  selected  at 
random,  say  2 99,  and  the  first  k group  is  sampled  at  the  299t.h  position. 

This  same  element  is  now  sampled  in  the  remaining  k segments,  i.e.,  numbers 
1099,  1899,  26 99>  and  3U99  on  the  original  numbering  scheme.  The  five 
samples  constitute  a set  of  systematic  samples  along  the  profile./ 

In  profile  D stratified  systematic  sampling  is  used.  When  the  strata 
ere  of  equal  length  it  is  only  necessary  to  select  a number  at  random 
within  the  range  of  sampling  positions,  and  to  repeat  the  s ampling  at  the 
same  point  in  each  stratum.  In  the  present  case  the  strata  are  of  unequal 
length,  and  it  was  decided  to  select  a number  between  00  and  99  to  represent 
a proportional  distance  in  the  stratum.  The  random  number  table  yielded 
3h>  which  was  used  as  0.3U  to  locate  a sample  at  that  proportionate 
distance  r**jm  the  initial  point  in  each  stratum. 

In  considering  these  four  alternative  methods  for  sampling  along  a 
profile,  the  question  may  be  raised  whether  any  one  of  the  methods  is 
superior  to  the  others.  If  the  beach  population  along  the  profile  is 
homogeneous,  simple  random  sampling  provides  a satisfactory  estimate  of 
the  beach  population  parameters.  If  the  beach  population  varies  markedly 
along  the  profile,  simple  random  samples  may  in  some  instances  give  a 
poor  estimate,  because  the  random  numbers  may  all  fall  into  a single  natural 
beach  zone  or  stratum.  Inasmuch  as  experience  amply  demonstrates  the 
occurrence  of  these  natural  zones  on  the  beach,  stratified  sampling  seems 
to  be  preferable.  Whether  samples  taken  at  random  in  each  stratum  are 
preferable  to  samples  systematized  in  each  stratum  probably  cannot  be 
answered  without  some  knowledge  of  the  population  gradients  involved. 

Cochran  (1953)  discusses  some  of  the  factors  involved  in  populations 
with  gradients,  with  cyclical  changes,  etc. 

When  a beach  area,  as  opposed  to  a profile  line,  is  considered,  the 
several  types  of  sampling  become  applications  of  plane  sampling,  also 
discussed  by  Cochran  (1953,  p.  183)  and  developed  in  greater  detail  by 
Quenouille  (19U9).  Suppose  an  area  100  x 100  feet  on  the  backshore  is  to 
be  studied  for  its  characteristics,  and  it  is  decided  te  take  16  samples 
from  the  area.  If  each  sample  is  to  occupy  a 3-inch  circle,  there  are 
some  160,000  possible  sampling  positions  available. 

The  100-foot  edges  of  the  large  square  are  numbered  from  0C1  to  IiCO 
in  sequence,  suggested  in  Figure  UA.  Sixteen  pairs  of  3-digit  numbers 
in  this  range  (omitting  duplicate  pairs)  are  drawn  from  a random  number 
table,  and  the  samples  are  taken  from  the  corresponding  positions,  using 
the  paired  numbers  as  x-  and  y-  coordinates.  This  yields  a set  of  16 
random  samples  from  the  beach  area. 

In  Figure  iiB  the  100-foot  square  is  divided  into  16  smaller  squares 
by  driving  stakes  at  25-foot  intervals  over  the  larger  square.  Each  of 
the  25-foot  squares  is  numbered  from  00  to  99  along  the  edges.  By 


drawing  16  pairs  of  numbers  in  this  range  and  collecting  samples  at  the 
corresponding  points,  a set  of  stfatified  random  samples  is  obtained. 

In  this  instance  the  subdivision  of  the  100-foot  square  into  smaller 
squares  is  a process  of  stratification  in  two  dimensions. 

In  Figure  l£  a sample  is  randomized  in  the  upper  left  25-foot 
square  by  drawing  a pair  of  random  numbers  between  00  and  99.  This 
same  position  is  now  sampled  in  the  remaining  15  squares.  The  16  samples 
constitute  a set  of  aligned  systematic  samples,  corresponding  to  stratified 
systematic  samples  along  a line.  When  samples  are  taken  directly  at  the 
stakes,  which  has  been  a common  procedure,  the  samples  have  the  same 
kind  of  pattern,  and  the  rnr..'omi*ing  element  is  preserved  by  randomizing 
the  position  on  the  backshore  of  one  corner  of  the  major  ino-fcct  square. 

Another  form  of  systematization,  yielding  unaligned  systematic 
samples,  is  defined  by  Quenouille  (19U9)  and  also  discussed  by  Cochran 
(p.  183).  This  is  shown  in  Figure  Id).  The  position  of  a sample  in  the 
upper  3e  ft  25-foot  square  is  randomized  by  selecting  a pair  of  numbers 
between  00  and  99 » representing  the  vertical  and  horizontal  coordinates. 

This  is  followed  by  randomizing  the  horizontal  coordinates  in  the 
remaining  squares  of  the  upper  row,  and  randomizing  the  vertical  coordinates 
of  the  remaining  squares  in  the  first  column.  The  resulting  pattern  of 
sampling  points  is  systematic  in  the  sense  that  the  distance,  horizontal 
or  vertical  between  samples  in  successive  strata  is  the  same. 

Quenouille  considered  the  relative  efficiency  of  these  several  sampl- 
ing plans,  as  well  as  several  others.  In  general  it  may  be  said  that 
if  the  backshore  area  is  homogeneous,  ary  of  the  methods  will  give  fair 
estimates.  If  a gradient  is  present  in  the  population,  however,  it 
seems  probable  that  the  ur.  aligned  systematic  samples  uill  detect  it  most 
effectively.  This  may  be  seen  by  comparing  Figures  I4C  and  I4D.  in 
a gradient  from  left  to  right,  for  example,  would  be  picked  up  at  the 
same  four  points  in  each  row.  In  id),  on  the  other  hand,  the  different 
positions  of  the  samples  in  each  row  would  result  in  16  intersections  of 
the  gradient. 

The  preceding  discussion  gave  equal  weight  to  each  3-inch  circle  in  the 
sampled  area,  which  would  require  rather  accurate  location  with  transit 
and  tape.  For  practical  purposes  one  may  decide  to  randomize  to  the 
nearest  square  foot,  taking  the  .sample  from  the  edge  or  center  of  the 
designated  square.  For  most  beach  studies  it  is  probably  adequate  to 
prepare  a drawing  or  map  of  the  sampling  design  iu  ‘die  office,  randomize 
the  sample  positions  on  the  map,  and  collect  them  in  the  field  to  the 
nearest  conveniently  measurable  unit. 

In  some  detailed  beach  studies,  and  in  connection  with  some  analysis 
of  variance  designs,  it  is  desirable  to  have  more  than  one  sample  in 
a grid  cell.  Multiple  samples  are  obtained  by  repeated  application  of 
the  same  randomizing  element  for  succeeding  samples.  In  similar  fashion, 
it  may  be  desirable  to  collect  more  than  oae  sample  from  each  beach  zone 
along  a profile.  This  may  also  be  accomplished  by  successive  randomi- 
zation of  the  sample  positions.  An  alternative  is  to  take  a number  of 
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samples  from  e ach  zone  proportional  to  the  zone  width.  This  g ives  rise  to 
a "represen+ative  sample"  in  that  each  zone  (or  stratum)  is  represented 
according  to  its  size. 


Multilevel  s ampltng  has  not  been  mentioned  explicitly,  although  in 
some  respects  it  is  implied  in  the  compromises  with  randomization  made 
in  an  waiiicr  raph.  Thus,  adoption  of  1 square  foot  as  a unit, 

sampling  area,  with  a 3-inch  sample  taken  from  the  center,  is  a form  of 
multilevel  sampling.  This  could  be  made  explicit  by  selecting  the 
square  foot  at  random,  and  then  randomizing  the  position  of  the  sample 
within  the  square  foot.  Whether  these  added  refinements  are  desirable 
in  routine  engineering  studies  may  be  questioned,  although  it  is  ap- 
propriate that  the  question  be  examined  academically. 

In  the  more  usual  form  of  multilevel  sampling,  one  might  select 
several  100-foot  squares  along  the  backshore  and  divide  each  into  25-foot 
squares.  Two  or  more  of  the  25-foot  squares  may  be  picked  at  random 
from  each  100-foot  square,  and  within  each  selected  25-foot  square  two 
or  more  samples  may  be  taken  at  random.  This  would  produce  a set  of 
samples  in  a three -level  hierarchy,  and  from  the  data,  the  variance 
contributions  of  each  level  of  s arapling  may  be  computed.  Inasmuch  as 
the  method  may  be  used  to  test  the  effects  of  sample  spacing,  an  example 
will  he  given  in  a later  section. 


Cluster  Sampling  is  a process  of  collecting  a group  of  closely 
spaced  samples  at  each  major  sampling  point.  Thus,  if  multiple  samples 
are  to  be  collected  from  each  beach  zone  along  a profile,  a stake  could 
be  driven  in  the  beach  at  the  randomized  sampling  points  (as- in  Figure 
3B),  and  two  or  more  samples  oould  be  collected  around  each  stake.  One 
way  of  doing  this  is  to  draw  a circle  about  2 feet  in  diameter  around 
the  stake,  and  divide  the  circle  into  say  16  equal  segments.  By  drawing 
numbers  from  01  to  16  (omitting  duplicates)  from  a random  number  table, 
the  sarrples  could  be  taken  from  the  corresponding  circle  segments,  say  1 
foot  from  the  stake.  To  avoid  finite  population  corrections  (Cochran, 
1953)  the  number  of s egments  should  be  fairly  large  compared  to  the 
number  of  samples  in  a cluster. 


Cluster  sampling  as  described  here  involves  randomization  of  the 
main  sampling  point  and  of  individual  samples  in  the  cluster.  The 
rand  sample  itself  is  a type  of  cluster  sample  in  which  the  randcmizin 
process  is  applied  to  the  total  sample,  but  not  to  the  individual  sand 
grains.  This  type  of  sample  is  sometimes  referred  to  as  a "complex 
random  sample". 


The  many  implications  of  the  several  sampling  plans  treated  here  are 
not  further  discus Bed  in  this  report,  although  Cochran  (1953)  develops 
the  corresponding  statistical  theory  for  each  kind  of  sample.  In  part 
the  selection  of  the  plan  depends  upon  costs,  and  upon  whether  the  main 
objective  is  to  obtain  a population  estimate  or  to  stucfy  relationships 
within  the  populatibns  more  analytically.  Beach  sampling  includes  both 
aspects,  and  Cochran's  book  is  recommended  as  a basic  text  for  further 
analysis  of  the  sample  data. 
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ANALYSIS  OF  BEACH  SAMPLES 

The  preceding  section  showed  that  general  principles  of  sampling  may 
be  applied  tc  beaches  and  that  the  application  of  these  principles  pro- 
vices  a number  of  alternative  ways  for  laying  out  samples  along  profiles 
or  on  a grid.  The  question  of  the  best  choice  of  method  depends  in  part 
on  the  objectives  of  the  s tudy,  and  this  in  turn  depends  upon  the  use 
to  be  made  of  the  samples,  in  this  section  some  of  the  kinuB  of  analyses 
performed  on  beach  samples  are  discussed,  to  pave  the  way  for  more  detailed 
consideration  of  the  selection  of  specific  sampling  methods. 

Staples  from  beaches  are  analyzed  for  various  attributes.  The  most 
common  is  particle  size  analysis,  which  may  be  conducted  with  sieves, 
with  a settling  tube,  or  by  measuring  grains  under  the  microscope.  There 
is  considerable  discussion  of  the  relative  merits  of  these  methods  in 
the  literature,  inasmuch  as  each  involves  a different  basic  definition 
of  "size”.  For  practical  purposes  it  is  perhaps  sufficient  that  all 
size  data  used  in  any  one  study  be  obtained  by  the  same  process,  so  that 
the  mean  sizes  obtained  represent  the  same  operationally-defined  variable. 
Only  when  sets  of  data  involving  different  methods  are  compared  is  it 
necessary  to  show  that  one  method  is  equivalent  to  another.  Various 
statistical  tests -for  showing  such  equivalence  or  its  lack  are  available. 

In  addition  to  variations  in  methods  of  analysis,  there  is  also  con- 
siderable variation  in  methods  of  expressing  mean  particle  size.  Perhaps 
the  most  commonly  used  is  the  median  diameter  in  millimeters,  which  re- 
presents the  size  of  the  middlemost  grain  in  the  weight  frequency  dis- 
tribution. Earlier  practice  had  been  to  use  the  arithmetic  mean  diameter 
in  millimeters,  graphically  determined  as  the  "equivalent  grade*.  Re- 
cognition that  many  beach  deposits  approach  log  normalcy  in  their  size 
distribution  led  to  the  development  of  the  logarithmic  mean  (phi  mean) 
and  its  antilog,  the  geometric  mean  diameter  in  millimeters.  Use  of  the 
log  median  (phi  ITK/Uidn)  15  CU.3G  COiiuTiOu  with  some  workers. 

In  similar  manner,  the  degree  of  spread  of  the  grains  about  the  mean 
size,  is  expressed  in  various  ways.  Those  most  commonly  used  are  the 
grading  factor, a variant  of  the  mean  deviation;  the  sorting  coefficient 
baaed  on  quartile  ratio;  and  a log  standard  deviation  (phi  standard 
deviation)  based  on  the  log  transformed  data. 

The  relative  merits  of  the  methods  of  expressing  particle  size 
are  not  of  major  concern  here,  except  to  point  out  that  in  general  the 
more  suitable  measures  are  those  in  which  the  mean  values  tend  to  approach 
a normal  distribution,  and  in  which  the  measures  of  spread  tend  to 
approach  a Chi-square  distribution.  Statistical  analysis  of  the  sample 
data  is  facilitated  under  these  conditions.  Fortunately,  the  mean  grain 
size  in  all  methods  of  analysis  is  based  upon  a large  number  (or  weight) 
of  grains,  so  thau  most  of  the  average  values  used  probably  tend  toward 
symmetrical  distributions.  This  point  needs  some  verification;  perhaps 
more  important,  the  sample  distribution  of  the  various  sorting  measures 
needs  critical  study. 
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In  any  event,  particle  size  analysis  is  usually  expressed  as  an 
average  diameter,  a measure  of  spread,  and  in  some  instances  a measure 
of  skewness*  These  are  the  data  that  are  generally  available  for  further 
statistical  analysis. 

The  fact  that  each  beach  sample  provides  an  estimate  of  average 
diameter,  spread,  and  skewness,  raises  certain  complexities  in  further 
statistical  analysis  of  a group  of  samples.  It  was  mentioned  that  beach 
sand  samples  are  a kind  of  cluster  sample  because  each  sample  contains 
a large  number  of  grains.  Inasmuch  as  the  randomizing  process  is  applied 
to  the  cluster  and  not  to  the  individual  grains,  the  latter  are  not 
independent  in  a statistical  sense.  Moreover,  the  mean  diameter  and  other 
measures  from  each  sample  are  commonly  based  on  weight  percentage 
frequency,  whereas  further  analysis  is  based  on  number  frequency.  Finally, 
the  sample  average  is  commonly  expressed  as  a median  diameter,  whereas 
subsequent  analysis  involves  the  arithmetic  mean  of  the  medians.  The 
assistance  of  mathematical  statisticians  will  be  required  to  sort  out  the 
variances  contributed  by  each  of  these  procedures. 

In  addition  to  particle  size  analysis,  beach  samples  are  commonly 
studied  for  their  mineral  composition.  Most  beach  sands  are  composed 
predominantly  01  quartz  and  other  light-colored  minerals,  with  a small 
percentage  of  darker  minerals.  These  darker  minerals  have  a greater 
specific  gravity  and  are  separated  from  the  sand  with  heavy  liquids. 

The  heavy  minerals  are  important  in  identifying  the  source  rocks  and 
source  areas  of  the  s and.  Common  practice  is  to  count  several  hundred  to 
a thousand  grains,  noting  the  abundance  of  the  several  species  present. 

The  results  are  commonly  expressed  in  number  percentage  terms.  As  was 
pointed  out  previously,  the  distribution  of  minerals  in  sediments  tends 
to  be  binomial  for  the  more  abundant  ones  and  Poisson  for  the  rarer  ones. 

In  some  special  studies  the  shape  of  the  grains  is  also  studied. 

Shape  may  be  expressed  in  tepms  of  the  sphericity  (operationally  defined 
as  the  cube  root  of  the  ratio  of  particle  volume  to  volume  of  a sphere 
circumscribed  about  the  particle),  and  as  the  roundness  (operationally 
defined  as  the  ratio  of  average  radii  of  corners  and  edges  to  radius  of 
circle  inscribed  in  maximum  projection  plane  of  the  particle).  Various 
short  methods  for  estimating  these  properties  are  available. 

Particle  orientation  (compass  direction  of  longest  particle  axis 
in  space,  for  example)  is  not  extensively  studied  for  engineering 
purposes,  although  some  methods  are  available  for  pebbles  and  large  sand 
grains. 

Other  kinds  of  observations  which  may  bo  made  on  beach  deposits  include 
study  of  beach  bedding,  beach  firmness,  beach  slope,  and  such  special 
features  as  beach  cusps,  etc. 

From  a statistical  point  of  view  the  size  of  a sample  should  be  ad- 
justed to  the  beach  characteristics  being  studied.  In  practice,  however, 
most  beach  samples  are  of  such  size  that  abundant  grains  are  available 
for  any  kind  of  analysis.  In  special  studies  it  may  be  found  that  sample 
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spacing  is  influenced  by  the  data  of  main  interest.  Decisions  on  this 
point  must  in  part  await  additional  knowledge  of  the  population  gra- 
dients of  particle  size,  particle  shape,  mineral  content,  and  the  like. 


Perhaps  the  most  common  problems  that  arise  in  connection  with  the 
several  objectives  listed  in  the  introduction  are  concerned  with  estimates 
of  the  over-all  characteristics  of  the  beach  deposit;  with  evaluation  of 
changes  in  texture  or  composition  along  and  across  beaches;  or  with  com- 
parison of  deposits  on  the  same  beach  during  different  seasons.  I tiaoimich 
as  average  particle  3izs  is  one  of  the  sain  considerations  in  such  comparisons 
it  is  used  here  as  the  main  illustration,  although  the  general  approach 
is  also  applicable  to  questions  of  source  areas  for  beaches. 


oujorAru-oui!  ur  odu  ijunvu  L/A1A 


In  proper  perspective  beach  sampling  is  a means  toward  an  end.  The 
objective  of  the  study  is  to  learn  something  about  a beach,  but  the  kinds 
of  tests  that  can  be  made  on  the  data,  and  the  convenience  of  computation, 
are  influenced  by  the  number  and  kinds  of  samples  available. 

As  a foundation  for  the  closing  section  of  this  report  on  recommended 
sampling  procedures,  the  beach  problems  mentioned  at  the  close  of  the 
preceding  section  are  treated  here.  Several  standard  statistical  methods 
are  introduced  to  indicate  way3  in  which  the  problems  may  be  approached. 

A basic  problem  that  arises  in  beach  studies  is  the  evaluation  of 
seasonal  changes  on  beaches.  Suppose  samples  are  collected  on  the  same 
beach  throughout  a year  and  it  is  desired  to  test  whether  there  is  a signifi- 
cant difference  in  average  particle  size  over  the  seasons.  Handin  (1951) 
furnishes  tables  of  data  which  may  be  used  as  illustrations.  Several 
samples  were  collected  each  month  from  a number  of  beaches.  When  several 
such  groups  of  samples  are  available,  it  is  possible  to  apply  analysis  of 
variance  to  the  data,  and  to  test  all  twelve  months  simultaneously. 

Certain  assumptions  underlie  applications  of  analysis  of  variance, 
as  discussed  by  Eisenhart  (1947),  and  Dixon  and  Massey  (1951).  Some 
preliminary  tests  applicable  to  the  data  are  mentioned  later. 

In  its  simplest  form  analysis  of  variance  compares  variances  between 
and  within  groups  of  data.  If  the  ratio  of  the  botween-group  variance 
to  the  within-greup  variance  exceeds  a selected  critical  value,  the 
difference  between  the  gioups  may  be  considered  statistically  significant. 
Dixon  and  Massey  (1951)  give  an  introduction  to  analysis  of  variance,  and 
describe  the  method  of  the  following  example  on  p.  121. 

Table  2 shows  the  data  arranged  in  month  columns.  The  nunber  of  ob- 
servations per  month  ranges  from  three  to  fiv».  The  design  used  is  a single 
factor  basic  form,  and  the  form  of  the  analysis  is  shown  at  the  bottom  of 
Table  2. 

Although  the  methods  of  computation  for  this  model  are  given  in  detail 
in  Dixon  and  Massey,  they  are  included  here  for  readers  unfamiliar  with  the 
technique.  The  computations  are  encumbered  by  extra  operations  because  of 
unequal  numbers  of  items  in  each  monthly  group.  The  first  step  is  to 
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square  all  the  items  in  the  table*  (321:)^  ♦ (U23)^  ♦ (3l*l)^  ♦ ♦ 

(335)2  ♦ (31tl  )^  ♦ (U66)2  = 6,1:25,  258.  This  may  be  called  Si  . The 
second  step  is  to  square  each  group  total  and  divide  it  by  the  number 
of  items  in  the  group,  and  then  add  the  quotients:  (l8lO)*/5  + (1910) */h  ♦ 

♦ (llU2)2/3  s 6,338,  U09.  Call  this  S2.  The  last  step  in  squaring 

is  to  square  the  grand  total  and  divide  it  by  the  total  number  of  items 
in  the  table:  (l69hh)2/U6  - 6,2Ul,286.  Call  this  S3. 

TABLE  2 

ARITHMETIC  MEAN  DIAMETER  OF  BEACH  SAND  COLLECTED  FROM  SAME 
BEACH  DURING  BACH  MONTH  (Data  from  Hand  in,  1951, 


Station  56,  p. 

108. 

Values 

w — * 

expressed  as  microns) 

Mar. 

Apr, 

Motr 
•— V 

June 

July 

h n rr  ^ 

Sep. 

Oct. 

Vmt 

»»w  . • 

i/w  a 

w «M»t 

Feb 

321: 

505 

355 

31:8 

291 

las 

311: 

350 

380 

392 

271* 

335 

1:23 

525 

370 

388 

3U5 

303 

319 

288 

383 

lao 

290 

3ia 

3ia 

1:60 

1*30 

322 

311 

1:58 

296 

388 

1:76 

382 

31:7 

1:66 

375 

1:20 

361 

1:32 

232 

318 

30U 

381 

3U7  U26 


TOTALS 

1810  1910  1516 

1916  91:7 

1U1 

At4f 

1620 

1^01  Al  1 

xluu  7J-X 

111:2 

T = 

362  1:78  379 

383  316 

353 

312  333 

1:05 

395  301: 

381 

ANALYSIS 

OF  VARIANCE 

Source 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean  Square 

F 

Between  Months 
Within  months 

97123 

8681:9 

n 

3U 

8829 

2551: 

3.1»6)h» 

Total 

183972 

1*5 

The  sums  of  squares  (SS)  in  Table  2 are  obtained  as  follows:  The 
total  SS  is  found  by  subtracting  S3  from  S]_:  6,1^5,258  ~ 6,2Ul,286  = 

323,972.  The  between-grcup  SS  is  round  by  subtracting  S3  from  S^s 
6,338,1*09  - 6,21:1,286  ■ 97,123.  Finally,  the  within-group  SS  is  obtained 
by  difference,  183,972  - 97,123  = 86,81:9. 

The  degrees  of  freedom  associated  with  the  groups  is  one  less  than  the 
number  of  groups,  12-1  - 11.  The  total  degrees  of  freedom  are  one  less 
than  the  total  number  of  items,  1:6  - 1 - U$.  The  within-group  sum  of  squares 
is  the  difference  between  these:  1:5  - 11  ■ 3U. 

The  mean  squares  are  found  by  dividing  the  sums  of  squares  by  the 
corresponding  degrees  of  freedom  to  obtain  the  values  shown  in  Table  2. 
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The  computational  process  thus  separates  the  total  variability  into  two 
sources  of  variability.  The  statistical  test  is  to  see  whether  the 
between-group  mean  square  is  significantly  larger  than  the  within-group 
mean  square.  In  making  this  test  one  has  to  decide  beforehand  on  the 
level  of  significance  at  which  he  is  willing  to  judge  the  results.  In 
the  present  instance  the  statistical  hypothesis  is  set  up  that  there  is 
no  significant  difference  between  the  group  means  at  the  conventional 
significance  level.  Selection  of  the  5%  level  implies  that  one  is 
taking  the  risk  of  falsely  rejecting  the  hypothesis  once  in  20  times  on 
the  average.  The  risk  of  making  this  f alse  decision  may  be  reduced  by 
accepting  a 1%  significance  level,  but  in  that  case  there  is  a greater 
risk  of  not  rejecting  the  hypothesis  if  the  means  really  are  different. 


In  the  present  experiment  the  level  is  chosen.  The  ratio  F of  the 
tween -group  "to  tlis  within-group  SQiisrs  is  computed  end  found  to  be 

3 c/46#  Tbe  critiCeij.  Value  of  the  F ratio  is  looked  up  in  a table  of  tbe 
F distribution  (Di^.on  and  Massey,  p.  310)  for  11  and  degrees  cf 
freedom.  By  interpolation  it  is  found  to  be  about  2.08.  The  observed 
value  is  larger  than  this,  and  hence  is  statistically  significant  at  the 
level.  The  hypothesis  of  no  significant  differences  among  the  group 
means  is  accordingly  rejected. 


One  may  also  check  the  F value  for  the  same  degrees  of  freedom  at 
the  1$  level.  The  table  on  p.  312  of  Dixon  and  Massey  yields  the  value 
2.82.  approximately.  This  is  also  smaller  than  3-li6,  so  that  the  ratio 
is  significant  at  the  l£  level.  This  is  shown  in  Table  2 by  adding  two 
stars  after  the  figure. 


One  may  conclude  from  this  experiment  that  there  is  a highly  signif- 
icant difference  between  the  monthly  means  of  average  particle  size  on 
this  beach.  This  implies  a strong  seasonal  change,  but  it  does  not 
identify  the  specific  months  in  which  the  changes  occurred.  Methods  of 
contrasts  are  available  for  such  tests,  but  a simpler  method  in  the 
present  instance  is  to  combine  the  original  data  into  two  main  groups, 
say  from  April  to  September  and  from  October  to  March  inclusive,  and  use 
the  t-test  on  the  two  groups  (Dixon  and  Massey,  1951,  p.  102).  When 
this  test  is  run  on  these  data,  it  is  found  that  a significant  gross 
seasonal  difference  is  also  present. 


The  one  factor  analysis  of  variance  form  nay  be  used  in  any  study 
involving  k groups  of  n items  each.  These  may  represent  samples  collected 
in  different  months  on  the  same  beach;  different  profiles  sampled  at  the 
same  tine  along  a given  beach;  samples  from  different  beaches  taken  at 
the  same  time;  or  groups  of  samples  from  different  zones  on  a given 
beach. 


A somewhat  more  versatile  analysis  of  variance  design  of  considerable 
value  in  beach  studies  is  the  row-column  (two  factor)  form.  This  permits 
simultaneous  comparison  of  two  factors,  such  as  changes  along  and  across 
beaches;  different  parts  of  the  s ame  beach  sampled  at  different  times; 
different  beaches  sampled  at  different  times;  and  other  combinations. 


/ 
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Fcr  tr:is  design  sa-scles  nay  oe  crllec'ed  on  Frits , and  the  fc Hewing 
illustration,  taken  fre-n  an  earlier  paper  bjr  t he  v ri ter  (1553)  illustrates 
a study  of  population  r.orocgeneity  on  a beach  at  Wilmette,  Illinois.  The 
data  are  shewn  in  Table  3«  The  upper  three  rows  of  sanples  were  taken 
frert  the  backshore,  and  the  lower  row  was  taken  fren  the  foreshore. 


TABLE  3 


o n 
e 
« 

23 

C 
«s 
O 
u 


GRID  CF 

VARIATE 

Gw  IK  t rc. 

P.&R.TI  CL 

v ct7c.  DATA 

EXPRESSED 

AS  PHI 

JEASS  OF 

PARTICLE  SIZE 

! DISTRIBUTION 

ALONG 

BEACH 

Backsbore 

2.15 

2.06 

2.10 

1.88  1.63 

Hack shore 

2.10 

2.06 

1.66 

1.95  2.06 

Beck she re 

1.96 

1.90 

2.0? 

1 _R7  2 OC 

Foreshore 

0.2? 

0.52 

0.9? 

1.18  1.21 

AJ4A1T5IS  Cr  VARIANCE 


So  irce 

Sun  of 
Squares 

Degrees  of 
Freedom 

Mean  Square 

He  tween  columns 

0.0766 

h 

0.0192 

i *gfll 

Between  rows 

s.0529 

3 

1.661*3 

(across  beach) 

Residual 

0.763ii 

12 

0.0636 

Total 

5.3929 

19 

< 1 J£S 
26.  LS*-' 


The  problem  here  was  t c ccracare  aver  a re  particle  size  along  and 
across  the  beach  on  a sampling  grid  that  straddles  the  hem  line.  The 
statistical  hypothesis  to  be  tested!  is  that  there  are  no  colimn  effe  ts 
(along  beach)  and  no  row  effects  (across  teach)  at  the  5%  significance 
level.  The  design  and  computations  are  described  by  Dixon  -and  Jtassey 
(p.  127).  The  total  variability  is  separated  into  three  parts,  the  first 
associated  with  the  column  effect,  the  second  with  the  rew  effect,  and 
the  third  with  a "residual " variance  which  represents  the  a r- or  terms . 

The  two  F ratios  are  ootainea  try  dividing  the  column  ar.d  row  mean  squares 
respectively  'ey  the  residua  I mean  square. 

Tbe  analysis.  shews  that  there  is  a higniy  significant  row  effect, 
but  no  significant  col  urn  effect  (as  indicated  by  the  letteis  ME). 

This  means  in  essence  that  the  backshcre  population  may  be  considered 
separate  from  the  foreshore  population s whereas  in  both  these  populations 
there  is  no  significant  change  in  means  along  the  length  of  the  grid 
(about  200  feet). 


when  prominent  zone  changes  occur,  aa  in  this  example,  the  use  of 
the  gross  average  to  characterise  beach  particle  size  is  generally  less 
effective  than  the  use  of  two  means,  one  for  the  backshore  and  ore  for 
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the  foreshore.  In  the  original  paper  the  separation  of  the  populations 
was  affected  by  a second  analysis  of  variance. 

On  the  basis  of  the  second  test  it  was  inferred  that  the  three  upper 
rows  ccnstitiite  one  population  and  the  lowest  row  another.  Hence,  in 
estimating  population  means,  the  values  in  the  upper  three  rows  may  be 
separately  averaged  to  obtain  a phi  mean  of  1.99*  corresponding  to  a 
geometric  m°an  diameter  of  millimeter;  The  foreshore  nhi  mean . 

on  the  other  hand,  is  0.83,  corresponding  to  a geometric  mean  diameter  of 
0.?6  millimeter.  Thus,  the  foreshore  particle  population  is  more  than 
twice  as  coarse  as  the  backshore.  A gross  average  of  all  values  in  the 
table  is  1.70  for  the  phi  mean,  corresponding  to  0.31  millimeter  for  the 
geometric  mean.  This  value  is  not  closely  representative  of  any  sample 
in  the  grid. 

A suggestion  from  this  example  is  that  if  estimates  are  to  be  made 
of  average  particle  size  on  given  parts  of  beaches,  the  sampling  plan 
should  involve  stratification.  By  including  in  each  group  only  samples 
from  a homogeneous  population,  the  mean  values  will  be  more  characteristic 
of  the  zone,  and  the  variance  among  the  means  will  generally  be  smaller. 

The  examples  presented  here  represent  several  basic  statistical  models 
for  comparing  groups  of  data.  More  advanced  analysis  of  variance  designs 
are  available  for  problems  of  wider  scope.  For  example,  use  of  a three- 
factor  analysis  of  variance  design  would  permit  simultaneous  evaluation 
of  changes  along  a beach,  across  a beach,  and  seasonally.  An  example  of 
this  type  of  design  is  given  by  Krumbein  (1953).  Dixon  and  Massey  (1951) 
is  highly  recommended  as  an  introductory  text,  both  for  analysis  of 
variance  and  for  definitions  of  terms  used  in  this  paper.  Heel  (392*7) 
provides  a somewhat  more  mathematical  introduction  to  statistics. 

The  analyses  of  variance  designs  included  here  are  used  mainly  for 
testing  hypotheses  about  population  means  on  a single  beach.  Yet  any  one 
beach  is  only  a sample  from  the  much  larger  population  of  all  beaches. 

If  interest  shifts  from  consideration  of  fixed  relations  among  means  to 
consideration  of  random  deviations  in  the  characteristics  of  any  one 
beach  from  the  mean  value  of  these  characteristics  in  the  much  larger 
population  of  all  beaches,  a model  appropriate  for  estimation  of  variance 
components  may  be  used.  An  xcellent  treatment  of  these  two  kinds  of 
analysis  of  variance  models  is  given  by  Eisenhart  (192*7). 

SAMPLING  PROCEDURES 


In  preceding  sections  no  question  was  raised  about  sample  size  ex- 
cept to  point,  out  that  usual  practice  affords  samples  of  adequate  size  for 
most  purposes.  Analysis  of  variance  may  be  used  to  study  this  problem 
more  critically,  and  elsewhere  the  writer  (1953,  1952*)  has  shown  that 
the  diameter  of  the  sample  is  of  less  importance  than  the  depth  of  penetra- 
tion for  both  sand  and  gravel  beaches.  The  experiments  cited  apply  to  a 
sand  beach  at  Wilmette  and  a gravel  beach  in  Evanston.  More  recently  the 
experiment  was  repeated  on  the  backshore  at  Ocean  Beach,  Maryland,  with 
the  same  results.  This  experiment  is  shown  in  Table  2*.  On  the  basis  of 
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these  experiments  it  seems  relatively  safe  to  infer  that  a sand  sample 
from  about  1.5  to  3 inches  in  diameter  is  satisfactory  for  particle  size 
analysis. 


TABLE  U 


CAMPLE  -SIZE  STUDY  OF  BEACH  SAND  FROM  BACK SHORE, 
OCEAN  BEACH,  MARYLAND  (DATA  EXPRESSED  AS  PHI 
MEAN  OF  PARTICLE  SIZE  DISTRIBUTION.) 


TV>nth 

oV 

Diameter  of  Sample 

Sample 

1.5* 

3" 

6" 

12" 

3" 

2.0p 

2.06 

2.05 

2. 05 

6" 

1.97 

1.98 

1.97 

1.98 

9" 

1.87 

1.88 

1.86 

1.86 

12" 

1.69 

1.6? 

1.67 

1.66 

ANALYSIS 

OF  VARIANCE 

Sum  of 

Degrees  of 

Source 

Squares 

Freedom 

Mean  Square 

F 

Between  columns 

0,0006 

0.00020 

3,33  NS 

(sample  area) 

Between  rows 

.3170 

3 

.10570 

176.2 

(sample  depth) 

Residual 

.0005 

9 

Total 

0.3181 

IT 

Common  practice  is  to  emphasize  sampling  of  the  upper  few  inches  of 
the  beach,  but  there  seems  to  be  little  organized  information  on  the  re- 
lative differences  that  would  show  up  seasonally  or  along  a beach  as  a 
function  of  the  depth  of  sample  penetration. 

That  the  problem  of  sample  depth  is  important  is  recognized  by  most 
workers.  It  is  not  unusual,  to  have  segments  of  the  beach  up  to  several 
feet  thick  cut  out  or  filled  in  near  the  berm  during  storm  cycles.  If 
samples  in  one  zone  penetrate  the  surface  of  a newly-deposited  layer, 
whereas  in  another  sons  the  samples  penetrate  a scoured  layer  normally 
several  feet  beneath  the  surface,  there  may  be  a reel  question  whether  the 
samples  are  strictly  comparable  in  terms  of  operating  processes.  There 
is  no  simple  way  of  avoiding  such  sampling  complexities,  although  some 
evidence  on  the  occurrence  of  marked  changes  could  be  had  by  running 
elevations  on  each  profile  at  successive  times  of  sampling  to  determine 
whether  any  serious  scour  had  occurred. 

There  is  some  justification  in  the  practice  of  sampling  the  upper 
few  inches  despite  these  complexities.  From  an  engineering  point  of  view 
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an  important,  cons  ice  ration  is  to  detemir.s  the  beach  characteristics  ir» 
the  zone  of  active  hydraulic  forces,  and  if  scour  ana  fill  occur,  that 
very  fact  jugr  add  to  the  variability  of  the  beach  surface  through  tire. 

S torn -cycle  or  seasonal  variability  nay  be  important,  in.  groin  design, 
fcr  example,  in.  that  changes  in  average  particle  ciz-e  on  the  foreshore 
through  the  seasons  will  be  accompanied  by  changes  in  foreshore  slope. 

Sore  knowledge  of  such  changes  is  valuable  in  designing  grcin.s  adjusted 
to  the  expected  range  in  foreshore  slope,  entirely  aside  from  the  specific 
processes  vhi/*h  way  produce  the  changes* 

Space  limitations  prevent  complete  discussion  of  sampling  devices, 
although  the1  subject  is  important  in  its-  influence  on  the  quality  end 
■quantity  of  s anple  obtained.  Samples  on  the  exposed  beach  na  be  collected 
by  digging  a hole  tc  a given  depth,  mixing  the  material  extracted,  and 
taking  a pint  or  quart  sartrle  by  seme  suitable  sub  sampling  7.ethod.  Per- 
haps more  cotcx-n  is  the  practice  of  rotating  a pint  or  half -pint  ice 
cream  carter  into  the  sand..  In  sene  instances  the  top  half- inch  cf 
surface  sand  is  first  scraped  off.  normally  s single  sample  is  taken  at 
each  sampling  point. 

Underwater  s an, piers  appear  to  be  of  three  raain  kinds.  The  sampler 
collects  a single  s angle  either  by  penetration  into  the  bottom,  by 
scraping  nr-  a surficial  layer,  or  by  "grabbing”  a sample  in  a snapping 
motion.  The  relative  efficiency  of  samplers  has  bees  locked  into,  and 
si era  could  be  done  along  this  lone  with  designed  experiments.  For 
present  purposes  it  will  be  assigned  tiiat  a sampler  is  used  which  provides 
samples  of  about  the  same  size  and  depth  of  penetration  as  are  collected 
on  the  above -water  part  cf  the  teach. 

Meat  underwater  samples  are  collected  along  traverses  in  terjrss  of 
water  depth.  The  samples  may  be  collected  at  increasing  depth  increments 
of  two  feet,  for  example,  tc  a total  depth  of  33  or  more  feet.  An 
alternative  method,  apparently  not  much  used  by  engineers,  is  to  space 
the  underwater  samples  according  tc  distance  from  the  near  tide  line, 
regardless  of  the  local  water  depth. 

If  comparison  of  underwater  samples  are  to  be  made  on  the  basis  of 
equivalent  water  depth,  scae centre!  of  this  factor  is  desirable,  Cn  the 
other  Hand,  if  it  is  true  that  the  nearshore  bottom  displays  cones 
similar  to  back  shore,  foreshore,  etc.,  it  may  be  that  sampling  within  these 
zones  will  provide  equally  useful  data.  The  writer  is  not  aware  that 
"popujatioa  zones'11  or.  the  nearshore  betto**  have  been  studied  as  such, 
although  daw  are  available  frum  aaps  that  shew  the  nearshore  distribution 
of  average  particle  sizes,  degrees  of  sorting,  etc.  Inman  (1953]  provides 
a number  of  such  maps,  and  recent  analysis  of  his  data  by  Rillar  (1952*) 
shows  that  certain  areas  can  te  distinguished  on  the  basis  of  multi- 
variate analysis  of  near,  size,  degree  of  sorting,  and  skewness. 

The  writer  has  observed  that  on  seme  Lake  Michigan  beaches  there  is 
a marked  coarsening  of  material  from  near  the  mean  water  line  to  the 
plunge  zone.  Beyond  this  the  bottom  is  score  tires  homogeneous  tc  moderate 


4 


f 


H 

i 


23 


depths.  Mats  by  Fisher  (1954)  an  cw  rough  h ne-ari  ty  of  underwater  cooes, 
3caetiz.es  at  an  angle  to  the  shore  line.  No  abreak"  as  conspicuous  as 
the  bera  for  separating  the  underwater  populations  is  present. 

The  problem  of  sample  spacing  arises  frequently  in  beach  sampling. 

Hop  far  apart  should  samples  be  spaced  along  a profile?  nos  far  apart 
should  the  profiles  themselves  be  spaced  along  the  beach?  Multilevel 
sane  ling.  renti cried  earlier,  provides  one  way  of  attacking  the  phoblen . 

Critical  data  or*,  beach  particle  site  are  not  available  as  an  ill- 
ustration, but  some  beach  firsness  data  obtained  at  Ocean  Beach,  Maryland 
for  another  purpose  in  1953  may  be.  used,  k portion  of  the  back  shore  is 
arranged  into  six  100-foot  squares  as  shewn  in  Figure  5 * Savh  of  these  is 
divided  into  nine  33-foot  squares.  Two  of  the  33-foot  squares  were  selected 
at  rand oas  in  the  100-foot  squares,  and  within  the  selected  33-fcot  squares 
two  rand  adzed  penetrometer  readings  were  sad  a.  The  two  values  in  each 
randomized  33-focrt  square  are  indicated  in  Figure  5.  Readings  were  made 
with  s dial  penetrosster  at-  6— inch  depth. 

If  the  experiment  had  been  designed  for  this  multilevel  study,  tha  10Q- 
f oot  squares  would  bs  spaced  randomly  over  the  baekshors.  The  number  of 
scalier  squares  within  the  103-foot  squares  would  be  increased  to  avoid 
necessity  for  finite  population  corrections.  The  example  will  be  treated 
as  though  it  had  been  properly  designed,  in  order  to  shew  tha  kind  of  data, 
that  nay  be  obtained1  from  such  a study. 

Preliminary  .analysis  of  the  large  squares  as  six  groups  of  four  items 
eaen  in  a single  factor  design  similar  to  Table  2 indicated  that  the  back- 
shore  in  the  six  large  squares  is  homogeneous  in  beach  firmness. 

For  present  purposes  tha  variance  ccatpnnsnts  of  the  three  sampling 
levels  are  ic^ortant.  These  were  found  by  use  of  a components  of  var- 
iance model  for  suit! level  sampling  (Cochran,  1953,  chapter  10;  Anderson 
and  Bancroft,  1952,  p.  325;  Olsen  and  Potter,  1954).  The  values  were 
obtained  as  follows: 


Source  Variance  Ccaqxraent 

Between  200- foot  squares  44.3 

Between  33-fact  squares  within  100-foot  squares  43.9 

Between  samples  within  33-fcct  squares  226.0 

The  largest  variability  occurs  between  individual  sanples,  whereas 
the  variance  c or. tri cations  of  the  33-foot  .and  103 -foot  squares  are  roughly 
the  same.  This  o.der  need  not  be  general,  so  that  the  f ulowlog  inter- 
pretations are  specific  instances  rather  than  generalizations.  Finite 
population  corrections  are  ignored,  as  stated  earlier. 

In  evaluating  the  penetrometer  data  in  terms  of  the  relative  effects 
of  saaple  spacing,  it  will  in  general  be  true  that  the  samples  within  the 


24 


/ 


25 


FIGURE  5 - BEACH  FIRMNESS  SAMPLING , OCEAN  CITY  , MD.,  1953 


33-foot  square?  will  be  closer  together  than  the  squares  themselves 
-re  in  the  largest  units.  In  the  present  instance  one  nay  say  that 
the  variability  between  the  closely  spaced  samples  within  a 33-foct 
square  is  greater  than  that  between  33"foct  squares  or  between  lOC-iwt. 
squares.  If  one  wishes  to  obtain  the  best  estimate  of  average  penetra- 
bility in  the  backshore,  the  question  nay  be  raised  whether  it  is  better 
to  take  nore  readings  close  together,  or  acre  readirgs  spread  over  the 
larger  sampling  units. 

It  can  be  shown  (Anderson  and.  Bancroft,  1952,  p.  326}  that  if  the 
number  of  100-foot  squares  is  a,  the  number  of  selected  33-foot  squares 
within  one  100- foot  square  is  and  if  the  amber  of  samples  vl'Uiin 
each  33-fcot  square  is  c,  then  the  variance  of  the  grand  mean,  V(x),  is 

V(T)  = (iiU.3/a)  ♦ (a3.9/ab)  ♦ (226.0/abc) 

where  the  Kaisers  in  the  numerators  are  taken  iron  the  preceding  com- 
ponents of  variance  table.  In  the  present  experiment,  a r 6,  b = 2, 
c s 2.  Hence, 

?(X)  2 7.33  * 3.36  * 9.1:2  = 20.1*6. 


The  variance  of  the  re  an  is  a measure  of  the  variability  to  which 
the  sample  mean  is  subject.  By  using  the  figure  20. 1*6,  it  is  possible 
to  see  hew  -ii derivative  sampling  plans  would  decrease  this  value-.  If 
the  number  cf  samples  is  to  be  doubled,  this  could  be  done  by  applying 
the  53*  plan  to  twice  as  many  100- foot  squares  on  the  backshore,  or 
by  taking  twice  as  many  readirgs  in  the  cresent  set  cf  33“ foot  squares. 
In  the  first  alternative,  a - 12,  b - 2 ad  c * 2. 
alternative  a ; 6,  b r 2,  and  c r it< 


By  inserting  these  values  in  the  expression 
alternative  yields: 


for  7(f: 


} , the  first 


7(f)  = (i»L. 3A2)  ♦ (ii3.9/2a)  ♦ (226.0/16)  r 10.23 


The  second  alternative  yields: 


V(X)  r ( isir . 3/6}  ♦ ( u3 • 9/12 ) * (226.0 /aS)  = 1>. 7; 


It  is  evident  from  tnese  values  that  Uwuhlli^  the  rusher  cf  largest 
sampling  units  will  cut  the  variance  of  the  wean  in  ha^f,  whereas  by 
doubling  the  number  of  lowest-level  sampling,  units  the  variance  is  only 
decreased  about  25  oer  cent. 


If  generalisations  are  permitted  frost  this  experiment,  it  would 
appear  that  a better  estimate  cf  the  rear,  arises  fro-  taking  additional 
samples  over  tie  larger  units  than  by  increasing  the  number  of  closely 
spaced  samples  in  the  smaller  units.  This  conclusion  was  alsc  readied 
by  Potter  and  Olson  (l95d)  in  their  study  of  cross-bedding  of  Pennsylvanian 
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sandstones. 


Although  this  example  freer:  multilevel  sampling  does  not  completely 
answer  the  question  of  sample  spacing*  it  dees  indicate  that  statistical 
methods  are  available  for  attacking  these  problems.  Sene  informal  ex- 
periments by  the  w riter,  based  largely  cn  student  problems,  suggested 
that  nc  significant  differences  in  average  particle  sire  arose  fro t.  grid 
sampling  on  15-foot  cusd  jb-foct  spacings*  These  distances  may  re  of 
the  same-  order  of  magnitude  in  a sampling  sense,  however,  and  it  may  be 
that  greater  differences  in  relative  spacing  would  give  acre  conclusive 
results. 
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The  preceding  sections  of  this  report  indicate  sene  cf  the  com- 
plexities involved  in  beach  sampling.  It  dees  not  appear  possible  at 
present  to  propose  a single  unified  sampling  plan  vhich  can  take  account 
cf  all  these  complicating  factors.  It  is  possible,  however,  to  suggest 
certain  procedures  for  beach  sampling  vhich  increase  the  reliability  of 
the  data  and  gre s bly  simplify  their  statistical  treatment.  The  following 
reooisne r. d a t ic ns , qualified  where  necessary , represent  the  writer's 
present  Judgment  on  problems  of  beach  sampling: 

(l)  Sample  Sine  and  Depth.  Samples  should  be  of  the  same  rise 
(vol  Ur.  or  weight  as  far  as  possible  for  any  given  study.  This  eliminates 
the  necessity  of  weighting  means  and  permits  direct  comparison  of  sets 
cf  data.  Sample  size  involves  two  aspects:  cross-sectional  area  of  the 
sample  and  depth  of  penetration. 

a.  The  studies  of  sample  size  given  earlier  indicate  that  for 
sand  samples  a circular  area  between  about  1.5  and  u inches  in  diameter 
is  satisfactory.  A variety  of  cylindrical  cardboard  cartons  in  this 
range  is  available.  Ice  ere  an  cartons  appear  to  be  particularly  suitable. 
For  gravel  the  sample  area  probably  needs  enlarger, era t if  depth  of  pene- 
tration is  to  remain  the  sane  for  both  sand  and  gravel.  On  gravel 
beaches  the  samples  nay  be  taken  to  greater  depths  to  compensate  for  area; 
judgment  rather  than  fixed  rules  is  indicated  at  this  stage. 

b»  For  most  sasd  beaches  a sarsple  about  2 or  3 inches  in  depth 
sse?ss  to  be  satisfactory  if  surficial  features  are  to  be  emphasized.  If 
the  upper  half -inch  surface  layer  Is  first  scraped  off,  this  should  be 
done  for  all  samples.  In  special  studies,  where  individual  sedimentation 
units  are  to  be  investigated,  the  methods  discussed  by  Otto  (1938)  are 
applicable.  For  studies  involving  a greater  depth  of  penetra*ir>n  none 
standard  methods  cf  extracting  the  total  sample  and  cf  splitting  it  tc 
size  should  be  adopted.  One  way  is  to  dig  cut  a cylinder  of  sand,  say 
about  6 inches  in  di  are  ter  and  a foot  deep,  spread  the  sample  on  a 
tarpaulin',  nix  thoroughly,  and  quarter  successively  until  a field  sample 
of  about  1 pint  is  obtained.  This  rer — servts  a form  of  s desampling, 
and  for  fundamental  studies  the  stati  cal  effects  of  such  subsaaapling 
should  be  evaluated. 


It  Fjr-r  be  pointed  cut  that  current  5 asrvpl ing  practices  yield  very 
large  samples  in  a statistical  sense.  One  cubic  centimeter  of  beach 
sand  has  about  60,000  grains,  so  that  the  -usual  half-pint  carter,  con- 
tains about  3 cr  k million  sand  grains,  and  weighs  about  50C  grams. 

It  is  likely  that  very  auch  smaller  sasples  could  be  used,  especially 
for  microscopic  methods  of  particle  size  analysis. 

(2)  Ranctosisatica.  A randomizing  element  should  be  made  ar,  ex- 
plicit part  of  all  sampling  procedures.  Tables  of  random  numbers  are 
■generally  nest  satisfactory  for  this  purpose,  although  coin  tossing  mar 
be  used  to  decide  between  alternatives.  As  was  pointed  out  the  randomis- 
ing eleaent  eliminates  personal  bias  and  provides  a basis  for  applications 
of  sampling  theory.  For  most  studies  it  is  possible  to  prepare  the 
randomized  locations  in  the  office,  so  that  the  actual  field  sampling 
proceeds  according  to  plan.  Randomization  of  underwater  samples  is  pro- 
bably most  conveniently  arranged  on  a systematic  plan  along  profiles, 

to  avoid  complicated  Maneuvering  of  the  ship. 

(3)  Sanple  Stratification.  For  rest  beard,  studies  a plan  of 
stratified  sampling  is  probably  preferable  to  simple  random  sampling. 

RfV! crh.  lOTl  t**?  ?CC-?F5p2. 1 <1  £?y  vsTi^iig  qttto- t»  -fHg  2nd! 

noting  the  presence  of  bents,  of  firm  and  soft  areas,  of  gravel  patches 
or  other  size  irregularities,  of  magnetite  or  other  placer  streaks,  and 
so  on.  These  observations  may  be  used  for  dividing  the  beach  into 
sampling  zones.  An  alternative  is  to  adopt  standard  procedures  based  on 
a bsck^tore  zone , upper  and  lower  foreshore  zones,  and  up  to  several 
nearshore  bottom  zones.  Each  sene  is  then  sampled  separately. 


Underwater  zones  are  sore  difficult  to  distinguish  without  some 
preliminary  sampling,  so  that  the  usual  procedure*  of  ~ asupling  by  »ater 
depth,  or  a modification  based  on  systeruatic  samples  outward  from  mean 
tide  level,  sees  adequate  in  the  present  state  of  knowledge. 


(ii)  Sample  Spacing  and  ^uafcer  :f  Semples.  Decisions  cn  these  two 
aspects  cf  beach  sampling  are  aneng  the  most  important  in  the  study. 
Considerations  of  cost,  plus  the  fact  that  data  frcra  beach  samples  con- 
stitute only  a small  part  of  the  total  observations  made  on  beaches 
for  engineering  studies,  require  that  as  few  as  possible  samples  be 
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mmber  cf  safeguards  he  wishes  tc  include  in  his  sampling  plan, 
following  considerations  w»v  he  taken  into  — r-.-mt  * 


The 


a.  Particle  size  and  mineral  data  (as  well  as  such  measurements 
as  proper ticn  of  acid  solubles)  are  based  on  large  samples  in  a 
statistical  sense,  and  hence  each  sample  value  has  a reasonably  high 
degree  cf  reliability,  assuming  a standard  method  of  analysis.  In  one 
sense,  therefore,  a single  random,  sample  from  each  beach  zone  along  a 
profile,  (or  a set  of  systematic  stratified  samples  ) nay  afford  in 
estiaBte  cf  the  sene  characteristics  adequate  for  many  studies. 


b.  Although  each  particle  size  analysis  provides  .an  estimate 
of  average  particle  size,  of  spread  (sorting; , and  of  skewness,  these 
values  dc  vary  within  the  zones,  and  a better  estimate  is  obtained  by 
taking  snore  than  a single  sample  in  each.  It  is  the  writer’s  opinion 
that  where  tine  and  cost  permit,  a minimum  cf  four  samples  iron  each 
beach  zone  nay  be  collected  in  order  to  have  a basis  for  estimating  the 
variability  cf  means  within  the  zones,  horecver,  if  cosipar Isons  are  to 
be  made  between  different  zones  along  a profile,  or  between  the  sane 
zones  cn  different  profiles,  sets  of  four  samples  permit-  convenient 
applications  of  analysis  of  variance. 


c.  Sore  beach  observations  represent  single  values  rather  the  a 
averages.  Such  are  penetrometer  readings,  for  example,  in  contrast  to 
particle  size  analysis.  In  general  these  single  values  will  fluctuate 
mere  widely  than  average  values.  If  pe  net  rose  ter  readings  are  to  be 
taken  along  a profile,  it  may  be  advantageous  to  take  several  readings 
at  random  in  the  vicinity  of  each  sand  sample  rather  than  a single  read- 
IThi-S  vi X 1 prcTidf  s.  Ttc  ca-wyio . Anoiiii1? o t s 
is  that  they  tend  tc  be  normally  distributed  even  when  the  single  ob- 
servations form,  a skewed  distribution.  Beach  characteristics  listed  in 
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based  on  a reasonable  number  of  observations  are  used  instead  cf  angle 
observations. 


d.  Applications  of  subsampling  (multilevel  sampling)  to  beach 
problems  deserve  much  more  attention  than  teas  been  accorded  then. 

Estimates  of  gross  beach  characteristics  for  extensive  areas  on  the  beach 
can  be  evaluated  core  thoroughly  when  the  variabi?  It?  Introduced  by  each 
level  of  sampling  is  included  in  the  study.  It  ---av  e found  in  sore 
instances  that  considerable  work  can  he  saved  fo~  a riven degree  of  re- 
liability by  using  Efultilevel  sampling  methods  instead  of  detailed 
sampling  at  the  lowest  level. 

e.  The  spacing  of  samples  along  a profile  has  been  touched 
upon.  If  the  zones  are  sampled  separately,  thr  mirier  cf  samples  from, 
each  could  be  kept  constant  regardless  of  zone  .loth.  As  an  alternative, 
the  total  number  of  samples  along  the  profile  could  be  allocated 
according  to  the  relative  widths  of  the  beach  zones*  Each  of  these 
methods  carries  certain  statistical,  implications  in  variance  computations, 
and  the  choice  of  method  may  depend  upon  whether  ar.  average  of  the  entire 
profile  .1-5  to  be  estinisted,  or  whether  subsequent  analysis  :is  carried 

odLCftig,  ZCiiC  • 

A related  question  is  the  spacing  of  profiles  along  the  beach.  The 
known  lesser  variability  along  than  across  beaches  suggests  that  profiles 
may  be  spaced  further  apart  along  the  beach  than  individual  samples  are 
spaced  along  a profile.  Perhaps  some  general  rule,  such  as  a spacing 
equivalent  tc  three  times  the  average  profile  length  may  be  applicable. 

In  order  to  have  sufficient  data  for  reliable  estimates  cf  beach  properties, 
s-ome  jwnlEwm  number  of  profiles,  such  as  six,  should  be  taken.  For  mere 
detailed  study  of  local  populations  the  principle  of  plane  sampling  dis- 
cussed earlier  may  be  applied. 


29 


f.  The  number  of  namples  necessary  to  test  hypotheses  concerning 
differences  among  mean  values  depends  upon  the  magnitude  of  the  differences 
that  it  is  desired  to  detect.  The  variance  of  the  sampling  distribution 
of  th8  mean,  X,  is  O'2  /N,  where  <r2is  the  population  variance,  and  N is 
the  number  of  samples.  The  variance  can  be  reduced  either  by  increasing 
the  nudber  of  samples  or  by  decreasing  <r2  . The  latter  may  be  accomplished 
by  increasing  the  precision  of  the  methods  of  analysis.  In  practice  some 
combination  of  these  possibilities  may  be  used..  Dixon  and  Massey  (1951, 
p.  121,  219-220)  discuss  methods  that  may  be  used. 

(5)  A'?.alysi3  of  Samples.  The  same  methods  of  laboratory  analysis 
should  be  used  throughout  any  one  study.  This  assures  that  the  numbers 
used  in  the  statistical  analysis  have  the  same  operational  definition. 

That  is,  the  process  of  acquiring  the  numbers  i«  kept  constant,  so  that 
no  variability  is  introduced  into  the  data  beyond  that  inherent  in  any 
single  method  of  laboratory  analysis.  If  operator  judgment  is  required 
in  the  analysis,  a single  analyst  should  be  used,  or  provisions  made  far 
evaluating  operator  effects.  (See  Griffiths  and  Rosenfeld,  1954>  on 
this  point) * 

Statistical  analysis  is  facilitated  when  measurements  are  normally 
distributed.  This  may  be  checked  by  arranging  a cumulative  distribution 
as  in  particle  size  analysis,  and  plotting  the  points  on  arithmetic  pro- 
bability paper.  In  analysis  of  variance  the  data  may  depart  from 
normalcy,  although  the  variance  of  the  groups  should  be  independent  of  the 
group  means.  This  can  be  checked  by  a scatter  diagram.  If  sufficient  rep- 
lication is  available,  the  variance  of  the  groups  should  also  be  checked 
far  homogeneity  (Dixon  and  Massey,  1951,  p.  147).  In  some  instances  (as 
with  counted  data)  transformations  are  appropriate  (Bartlett,  1947,  and 
Kempthome,  1953,  p.  153)  • Table  1 may  be  checked  for  suggestions  on 
distributions. 

As  a general  rule  it  is  better  to  analyze  each  beach  sample  separately 
and  to  combine  the  data  algebraically,  than  to  mix  the  samples  before 
analysis.  This  separate  treatment  provides  data  for  estimates  of 
variability  between  samples,  which  is  lost  if  only  the  composite  data  are 
available. 

(6)  juiu^xju  .uxa.ov«^ovi  yja  a gxven  i «— 

present  only  one  of  a~large  number  of  possible  conditions  which  may  occur 
or.  beaches.  If  an  estimate  is  to  be  made  of  the  average  composition  of 
a beach  during  the  entire  year,  samples  will  be  required  during  the 
several  seasons  as  well  as  during  quiet  and  storaty  conditions.  It  is  not 
difficult  to  obtain  such  data  on  the  exposed  part  of  the  beach,  but  the 
collection  of  underwater  samples  may  not  be  feasible  during  severe  storms. 

Repetitive  sampling  may  be  designed  in  several  ways  to  obtain  a set 
of  samples  representative  of  the  range  of  seasonal  and  storm  cycle  con- 
ditions on  a beach.  Even  without  specific  seasonal  design,  sets  of  four 
samples  from  each  zone  collected  at  intervals  through  the  year  provide 
data  for  comparisons  by  analysis  of  variance  and  ether  techniques.  It 
is  probably  preferable  to  randomize  the  specific  sampling  positions  along 
the  profiles  separately  for  each  sampling  date. 
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Problems  of  weighting  the  data  may  be  involved  in  undesigned 
repetitive  sampling,  unless  some  provision  is  made  to  prorate  the 
number  of  sampling  intervals  according  to  seasons  and  storms.  For  ex- 
ample, if  1$  sets  of  samples  are  taken  on  quiet  days  and  only  one 
set  on  a stormy  day,  the  mean  value  may  not  be  characteristic  of  the 
average  beach  conditions.  Stratification  may  be  applied  to  the  problem 
to  overcome  some  of  these  difficulties.  Two  possibilities  appear  to 
apply: 


a.  The  year  may  be  stratified  into  months  or  seasons,  with 
samples  collected  systematically  or  at  random  in  each  time  stratum. 

Suppose  the  year  is  stratified  into  12  months,  with  two  sets  of  samples 
to  be  collected  from  each  zone  each  month.  The  s ampling  dates  in  each 
month  may  then  be  separately  randomized,  or  the  dates  may  h«  randomized 
for  the  first  month,  and  repeated  as  systematic  stratified  sampling  dates 
in  the  remaining  months.  An  alternative  here  is  to  stratify  according 

to  four  3-month  seasons,  with  six  randomized  or  systematized  sampling 
dates  in  each  season. 

b.  Instead  of  stratifying  the  year  according  to  months  or 
seasons,  the  total  year  may  be  divided  into' several  strata  representing 
quiet  intervals,  mild  storms,  severe  storms,  etc.,  according  to  data 
available  otr  the  relative  proportions  of  time  that  each  set  of  conditions 
may  be  expected  to  occur.  Then  by  distributing  the  2 U sets  of  samples 
over  these  strata  in  random  or  systematic  manner,  the  sample  data  will 

be  representative  of  physical  conditions  rather  than  simply  of  time. 

Stratification  of  the  year  into  time  periods  or  into  strata  defined 
by  occurrence  of  given  processes  involves  problems  of  cyclical  phenomena. 
Complete  randomization  within  the  time  strata  may  yield  more  representative 
samp}es  than  systematization  within  strata,  if  by  chance  the  systematic 
points  fall  at  corresponding  highs  or  lows  in  the  cycle.  Cochran  (1953) 
discusses  seme  of  the  factors  involved,  and  it  is  likely  that  collabora- 
tion of  mathematical  statisticians  will  be  required  to  solve  the  problem 
of  reliable  repetitive  sampling. 

CONCLUDING  REMARKS 

The  discussion  of  beach  samoling  in  this  report  is  given  from  a 
subject-matter  viewpoint  rather  than  from  a theoretical  statistical 
viewpoint.  The  recommended  sampling  procedures  are  the  writer's 
interpretation  of  sampling  principles  and  methods,  based  largely  on 
Cochran's  (1953)  classification,  and  on  a number  of  beach  sampling  experi- 
ments at  Nofcthvestem  University,  many  conducted  by  graduate  students. 

The  writer  has  no  misconceptions  about  the  tentative  nature  of  the 
sampling  methods  suK"|e.:ced  in  this  report.  The  p.cblem  of  beach  sampling 
requires  the  acti'ru  collaboration  of  mat:r?jn»-<  * cal  statisticians,  but  the 
statistician  canoe  of  little  help  if  poorly  collected  data  are  submitted 
to  him  and  he  is  what  can  be  done  with  them.  Sampling  design  is 

necessary  if  rigo  analysis  of  the  results  is  to  be  had.  It  is  the 
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writer's  hope  that  by  including  more  explicit  design  elements  in  future 
sampling  programs,  material  will  be  made  available  by  engineers  and 
geologists  for  further  analysis  by  statisticians.  From  the  results  of 
such  analysis,  modified  where  needed  by  practical  considerations,  valid 
sampling  programs  applicable  to  specific  beach  problems  can  be  developed. 
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